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ABSTRACT 

Several gravitationally lensed quasars are observed with anomalous magnifications in 
pairs of images that straddle a critical curve. Simple theoretical arguments suggest 
that the magnification of these images should be approximately equivalent, whereas 
one image is observed to be significantly demagnified. Microlensing provides a possible 
explanation for this discrepancy. There are two key parameters when modelling this 
effect. The first, the f raction of smooth matter i n the lens at the image positions, 
has been explored by ISchechter fc Wambsgansi (|2002l ). They have shown that the 
anomalous flux ratio observed in the lensed quasar MG 0414+0534 is a priori a factor 
of 5 more likely if the assumed smooth matter content in the lens model is increased 
from 0% to 93%. The second parameter, the size of the emission region, is explored 
in this paper, and shown to be more significant. We find that the broadening of the 
magnification probability distributions due to smooth matter content is washed out 
for source sizes that are predicted by standard models for quasars. We apply our 
model to the anomalous lensed quasar MG 0414+0534, and find a 95% upper limit 

of 2.62 X 10^^/iyQ^^^(M/MQ)^/^cm on the radius of the Tband emission region. The 
smooth matter percentage in the lens is unconstrained. 
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1 INTRODUCTION 

Gravitational microlensing, observed in some multiply- 
imaged quasars, oflters an excellent opportunity to constrain 
the size and structure of the regions responsible for emitting 
the radiation. Some lensed quasars exhibit anomalous flux 
ratios between pairs of images straddling a caustic that can- 
not be convincingly explained using simple macrolens mod- 
els. However microlensing models, particularly those that in- 
clude a significant smo oth matter component in the le ns, of- 
fer a possible solution (|Schechter fc Wambsgans3l2002l . here- 
after SW02). 

The quadruply imaged quasar MG 0414+0534 is one 
example of a lensed quasar displaying anomalous flux ratios. 
These anomalies are observed in ratios between images A2, 
located at a saddle poi nt in the time delay s urface , and A\, 
located at a minimum. ISchechter fc Moord (Il993l ) reported 
an A2/A\ I-band flux ratio of 0.45 ± 0.06 in observations 
taken on 2-4 November 1991. This result was supported by 
subsequent observations: A2 /A1 = 0.3 ±0.1 from CFHT in 
the I-band on 1 March 1992 (|Angonin-Willaime et "aLlll994l ') 
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and A2/A1 = .47 ± 0.01 from HST in the I-band on 8 
November 1994 (|Falco. Lehar fc Shapiroll 19971 ). 

However, lensing theory tells us that image magnifica- 
tion scales as the inverse of perpendicular distance from a 
critical curve l|Chang fc Refsdalll 19791 , blandford fc NaravanI 
| l986j ). Models that fit the observed image positions in MG 
0414+0534 place images A\ and A2 either side of such a crit- 
ical curve, and thus we would naivel y expect their magnifica- 
tion ratio, or flux ratio, to be ~ 1 (I Witt. Mao fc Schechtej 
1 19951 . hereaft:er WMS95). Indeed an 8 GHz radio flux ra- 
tio of ^2 Ml = 0.90 ± 0.02 was observed on 2 April 1990 
(|Katz fc HewittI Il993l '). This radio observation was much 
closer to the A2 /Ai ~ 1 ratio predicted using lensing models 
(WMS95). 

A mix of smooth and clumped matter distributions of- 
fers a potential explanation for the discrepancy between op- 
tical and radio flux ratios. Microlensing simulations that 
assume all matter in the lensing galaxy to be in compact 
objects yield a probability of 0.068 for a flux ratio lower 
than A2/A1 = 0.45 ± 0.06 (WMS95). However, the addi- 
tion of a smooth matter component in the lensing galaxy 
can significantly increase this probability up to values as 
large as 0.35, for sources with a characteristic radii much 
smaller than an Einstein Radius (SW02). The most com- 
pelling alternative explanation is millilensing by CDM sub- 



2 N. F. Bate, R. L. Webster and J. S. B. Wyithe 



structure jMetcalf fc Madaul l200ll . iMetcalf fc Zhaol |2002| . 
iDalal fc KochaneklbOoi T However millilensing should also 
eflFect the radio emission, for which an anomalous flux ratio 
is not observed in MG 0414+0534. 

Microlensing simulations have been used to place lim- 
its on the size of different emitting regions in several 
quasars. The majority of these efforts have focussed on 
the continuum region of Q2237+0305, and have relied 
on timescale arguments for high ma g nifica ti on events (eg, 
Warribsganss. Paczynski fc Schneider! Il990l . IWebster et al.l 



Table 1. Lensing paramaters 



199ll . IWvithe et al.ll200oi r These arguments depend upon 



the apparent transverse velocity of the source, and stellar 
proper motions within the lens, neither of which are well 
known. 

Upper limits on the size of emission regions can also be 
made using single observations of flux ratios between lensed 
images. These limits are free of the velocity dependency that 
arises in timescale arguments. For example WMS95 placed 
an 95% upper limit of 3 X 1O^'^(M/M0)^/^ cm on the I- 
band continuum emission regi on in MG 0414+0534 using 
this m ethod. On larger scales, IWavth. O'Dowd fc Websteil 
(|2005l ) placed an 80% upper limit on the broad emission 
line region in Q2237+0305 oi 2x10^'^ (M/Mq)^^^ cm. In both 
cases these limits are large enough that the emission regions 
cannot be considered as point sources in microlensing simu- 
lations. 

In this paper, we extend the investigation presented in 
SW02 to include the effects of varying source size on the 
magnification distributions of minimum and saddle point 
macroimages. We then use the anomalous observed flux ratio 
to place constraints on the size of the I-band emission region 
in MG 0414+0534, taking into account a varying smooth 
matter component in the lens. 

Our simulation method is discussed in Section 2, fol- 
lowed by a qualitative examination of magnification his- 
tograms in Section 3. In Section 4 we use our microlensing 
simulations to place upper limits on the size of the I-band 
emission region in MG 0414+0534. We discuss the results of 
this investigation in Section 5. 

Throughout this paper we use a cosmology with Ho = 
70fcms"Mfpc"\ fi™ = 0.3 and = 0.7. 



Image 


Type 


Ktot 


7 


IJ-tot 


MIO 


minimum 


0.475 


0.425 


10.5 


SIO 


saddle 


0.525 


0.575 


-9.5 


Ai 


minimum 


0.472 


0.488 


24.2 


Aa 


saddle 


0.485 


0.550 


-26.8 



Lensing parameters for generic low magnification minimum 
and saddle point ma croimages (MIO and SIO respectively, 
ISchechter fc Wambagansj 12003 ) , and for the im ages of interest 
in MG 0414+0534 (Ai and ^c. lWitt et"allll995l) . 

Magniflcation maps were generated with a resolution of 
2048 X 2048 pixels, covering an area of 247)0 x 24770, where 
770 is the Einstein Radius projected on to the source plane 
(3.75 x lO^'^ h~Q^^ {M/MQ^^^cm for MG 0414+0534). Each 
pixel therefore has a side length of O.OI770, which corresponds 
to 4.39 X 10"/ito ^^(WA^0)^''^cm for MG 0414+0534. 

The number of rays shot per unlensed pixel was chosen 
to give smooth probability distributions down to low mag- 
nifications. For the generic images 1000 rays per unlensed 
pixel were sufficient, whereas the higher magniflcation MG 
0414+0534 images required 4000 rays per unlensed pixel. For 
each model, 20 magniflcation maps were generated per im- 
age. This number of 24770 x 24770 maps was found to provide 
enough statistically independent data points to construct 
reasonable probability distributions. 

The magnification maps were then convolved with a 
Gaussian source intensity profile, with characteristic ra- 
dius varying from 0.05r;o (close to the source size used 
in SW02) to 2.007;o in steps of O.O5770. In physical units 
for MG 0414+0534, this corresponds to a characteristic 
source radius of 1.87 x IO^'^/ito ^^(M/M0)^''2cm to 7.50 x 
W^^h:^Q^^{M/MQy^^cm. Magnifications from the saddle 
point image were divided by those from the minimum image 
to construct flux ratios, or equivalently changes in magni- 
tude using Am — 2.5logio{iJ.2/ fJ,i)- Probability histograms 
were then derived from these changes in magnitude for each 
combination of smooth matter percentage and source radius. 



2 SIMULATIONS 

Microlensing simulation s were conducted using a 



rayshoot ing method (eg, Kayser, R efsdal fc Stabelll Il986l . 
IWambs ganss. Paczynski fc K atz 1990:). The key parameters 
for such simulations are the convergence utot and the 
shear 7 of the lens at the image locations. As discussed 
in SW02, the convergence of the lens can be split into 
two components - a continuously distributed component 
ACc, and a compact stellar component k,. We allowed the 
smooth matter percentage to vary from 0% to 99%. 

We conducted simulations for two sets of lensing pa- 
rameters (see Table[l|. The flrst (labelled MIO and SIO, fol- 
lowing SW02) represent generic minimum and saddle point 
images with total magniflcations /i ~ 10. These parameters 
were chosen both for their computational simplicity, and to 
allow comparison with SW02 results. The second set of pa- 
rameters, taken from WMS95, describe the Ai and A2 im- 
ages in MG 0414+0534. 



3 FLUX RATIO DISTRIBUTIONS 

Fig. [T] shows differential probability distributions of change 
in magnitude Am for 0% smooth matter content. Distri- 
butions are provided for the generic images MIO (mini- 
mum, left column) and SIO (saddle point, right column) and 
four characteristic source radii (0.05r;o, O.25r7o, O.5O770 and 
2.00r?o). 

For a source radius of O.O5770, the minimum and saddle 
point distributions are essentially identical. As source radius 
is increased both distributions become narrower, while pre- 
serving their overall shape. This is expected; fluctuations in 
source magniflcation due to microlensing are strongly de- 
pendent on the size of the region being le nsed, as has been 
noted many times in th e literature (eg, Chang fc Refsdal 
19791 IWambsganss et aLlll990l. IWvithe. Webster fc Turneii 



20od . iMortonson. Schechter fc Wambsgansj|2005l ). 

When smooth matter percentage is increased to 80% 
(Fig. [2}, the behaviour described in SW02 becomes ap- 
parent. For a characteristic source radius of 0.05r;o the 
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Figure 1. Differential probability histograms for Am for 0% 
smooth matter percentage. Histograms are displayed for images 
MIO (left) and SIO (right), and four characteristic source radii 
(top to bottom). The dashed line indicates the predicated macro- 
model magnification in magnitudes. 



Figure 2. Differential probability histograms for Am for 80% 
smooth matter percentage. Histograms are displayed for images 
MIO (left) and SIO (right), and four characteristic source radii 
(top to bottom). The dashed line indicates the predicated macro- 
model magnification in magnitudes. 



minimum and saddle point histograms are quite differ- 
ent. Both distributions are bifurcated, with a second 
peak towards Am — (the unlensed case). Separate 

J leaks correspond to different numbers of micro-images 
Granot. Schechter fc Wambsgans"3l2003l ). In addition, the 
saddle point distribution covers a much broader range in 
Am than the minimum distribution. 

Increasing the source radius in the 80% smooth matter 
case has a more extreme effect on the shape of the prob- 
ability distributions. For a characteristic radius of 0.25r;o, 
the second peak in both distributions has disappeared. A 
long tail towards Am — remains in the saddle point dis- 
tribution. This tail is suppressed by the time a character- 
istic source radius of O.5O770 is reached. As in the 0% case, 
increasing the source radius causes the distributions to nar- 
row. However, at source radius 2.00r7o the 80% smooth mat- 
ter probability distributions are narrower than in the 0% 
smooth matter case. 

We find similar behaviour in the 93% smooth matter 
case (Fig. [31). For a small source, minimum and saddle point 
distributions appear quite different. As the source size is in- 
creased, the distributions quickly become narrow. The fiat 
tail towards Am = in the O.O5770 saddle point distribu- 
tion is strongly suppressed by the time a source radius of 
O.5O770 is reached. The combination of a high smooth mat- 
ter component and a large source size make magnification 
fiuctuations due to microlensing very minor. 

These results differ from the predictions made in SW02. 
Although SW02 didn't present a detailed examination of 



the response of magnification histograms to increasing 
source size, they predict that doing so would compress 
the magnification histograms horizontally, while preserv- 
ing their shapes. On the other hand. Fig. [2] and Fig. [3] 
show quite clearly that increasing source radius compresses 
the probability distributions, and washes out the broad- 
ening of the distribution for negative changes in magni- 
tude. ICongdon. Keeton fc Osmer (|2007l ) observe the same 
behaviour for a different region of k — 7 parameter space. 
As a result, a large smooth matter percentage can only ex- 
plain anomalous flux ratios if the source size is small relative 
to the Einstein Radius. 



4 LIMITS ON SOURCE SIZE IN MG 0414+0534 

We repeated the above analysis for microlensing parameters 
appropriate for MG 0414+0534 and observed similar results. 
Fig. ^illustrates the effect of smooth matter and source size 
on the probability of observing anomalously low flux ratios 
in MG 0414+0534. For a small source radius, the Am cu- 
mulative probability distribution is considerably broader in 
cases where the smooth matter percentage is high. Given 
a small source, a flux ratio as low as the observed ratio 
{A2/Ai)obs = 0.45 has a probability of 0.07 for the 0% 
smooth matter case, and 0.31 for the 93% case. At a source 
radius of 0. 25770, however, the 0% and 93% smooth mat- 
ter distributions are virtually identical, and for larger radii 
the 93% smooth matter distribution is found to be narrower 
than the 0% smooth matter distribution. 



4 N. F. Bate, R. L. Webster and J. S. B. Wyithe 



source radius: 0.05?7q source radius: 0.2597q source radius: 0.50?7q 




Figure 4. Cumulative probability distributions for Am between images Ai and A2 in MG 0414+0534. Smooth matter percentages of 
0% (solid line) and 93% (dashed line) are displayed, for three characteristic source radii (left to right). The line representing the 0% 
smooth matter, 0.05t;o source radius distribution is provided in the second and third panels (dotted line) for comparison. 




Figure 3. Differential probability histograms for Am for 93% 
smooth matter percentage. Histograms are displayed for images 
MID (left) and SIO (right), and four characteristic source radii 
(top to bottom). The dashed line indicates the predicated macro- 
model magnification in magnitudes. 



We next compare the observed A2/A1 flux ratio for 
MG 0414+0534 with our model probability distributions 
and construct an a posteriori probability distribution for 
characteristic source radius and smooth matter percentage. 

Following WMS95 and SW02 we conduct our analysis 
using an observed flux ratio of Robs ~ (^2/^1)063 ~ 0.45 ± 
0.06 jSchechter fc Moordl 19931 ). By comparing this observed 
flux ratio with conditional probability distributions for the 
flux ratio, we constructed likelihoods for the observed ratio 
given varius source radii 77, in units of the Einstein Radius 770, 



and smooth matter percentages s = nc/ntot- Using Bayes' 
theorem, this likelihood L[Robs\s,Tf) was converted to an 
a posteriori differential probability distribution for smooth 
matter percentage and source radius as a fuction of Roba- 



dsd-q 



HP dP 

r / T-1 I \ przor i-tJ prior 

oc L{Robs\s, rj)- 



ds 



drj 



(1) 



These distributions are then marginalised over the observed 
distribution for flux ratio. 



d-'P 
dsdrj 



with a ■ 



dR 



d^P 



\ dsdr} 



2TTARob 



-(R — Rob 
■2^Rlb. 



(2) 



(3) 



where the error in the flux ratio was treated as a Gaussian 
with characteristic radius equal to the observational error 
of ±0.06. We used a logarithmic Bayesian prior for source 
radius, 

dP rior 1 

— oc — where i] ^ 2.OO770 
dri rj 

= otherwise (4) 
and a constant Bayesian prior for smooth matter percentage. 



dP 

i-tJ prior 

ds 



OC 1 where sC s < 0.99 
— otherwise 



(5) 



Source size is a quantity with units (in this case, units 
of Einstein Radius) , and is therefore assumed to have a prior 
probability that is constant per unit logarithm. The choice 
of a prior that is flat in the logarithm ensures that the ratio 
of prior probability for two values of source size does not 
depend on the units chosen. On the other hand, the value of 
the smooth matter fraction is a dimensionless quantity and 
is therefore assumed to have a uniform prior. 

Probability contours were then drawn through the re- 
sulting distribution and are plotted in Fig. \S\ This shows 
that our simulations do not constrain the smooth matter 
content of the lens. On the other hand an upper limit on 
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Figure 5. Probability distribution for characteristic source ra- 
dius, in units of the Einstein Radius rjo, and smooth matter per- 
centage. The contours show levels at 64%, 26% and 14% of the 
peak likelihood. These contours were constructed for an observed 
flux ratio of {A2 /Ai)ct,s = 0.45 ± 0.06 and a logarithmic source 
radius prior. 
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Figure 6. Cumulative probability that the I-band emission re- 
gion in MG 0414+0534 is smaller than a radius r\, given an ob- 
served flux ratio of {A2lA\)aba = 0.45 ± 0.06. Results are dis- 
played for two source radius priors: logarithmic (solid line) and 
uniform (dashed line). 



tical confidence of 95%. In physical units, this limit is 
2.62 X XQ^^h^^''^ (M IM^fl'^ cm. Our limit on the I-band 
emission region is smaller than the limit found in WMS95. 
The improvement is due to the inclusion of a variable smooth 
matter component in our lens. 

To test the sensitivity of our results to the assumed 
prior probability on source radius we have re-calculated the 
source radius constraints assuming a flat, rather than a log- 
arithmic, prior. These results are presented in Fig. [5] (dashed 
line). We find that the constraints are less stringent where 
the flat prior is used, with the data constraining the source 
radius to be smaller than about 1.3?7o. Thus, while the data 
do constrain the upper hmit on source radius, the precise 
constraints are sensitive to the prior chosen, indicating that 
the range of source radii that are consistent with the data 
remains considerable. 

We have assumed a fixed microlens mass of I.OMq 
throughout this analysis. Previous studies have found the 
mean microlens mass to be a good approximation to the full 
microlens mass fuct ion (eg, IWitt. Kayser fc Refsdall Il993l . 
iLewis fc Irwinll 19951 . IWvithe fc Turnedl2001^ . Under this as- 
sumption, a Salpeter mass function with a mass range of 
0.1 < M/Mq < 10 will reduce the physical size of our limit 
by a factor of ~ 2, or more if the lower mass bound is de- 
creased. 

However, recent studies have shown that the interaction 
between source size a nd microlens mass function is complex. 
ICongdon et"ai] ||200'<D found that including a mass spectrum 
for the microlenses leads to broader probability distributions 
for larger source sizes. The range of the mass function is 
important; its slope is not. This will lead to a relaxa t ion o f 
our upper limit on source radius. iLewis fc Gil-Merinol (|2006l ) 
demonstrated that small mass microlenses can mimic a 
smooth matter component if the source size is large enough. 

We have used a Gaussian source intensity pro- 
file throughout our analysis. This is reasonable as 
iMortonson et all (|2005l ) have shown that microlensing fiuc- 
tuations are insensitive to all properties of the source model 
except the radius. We note, however, that in a sheared mag- 
nificatio n map source elliptici ty and position angle are im- 
portant jCongdon et al.ll2007^ . 



the size of the I-band emission region exists for all smooth 
matter percentages. 

We can therefore place a limit on the size of the I-band 
emission region in the quasar. To do so, we marginalise the 
differential probability distribution over smooth matter per- 
centage: 



dr] J dsdri 

We then find the probability that the source is smaller than 
a particular radius, given the observed flux ratio: 

'7 

P«.)=/grfV (7) 



Fig. |6] shows the results of this analysis for an ob- 
served flux ratio of (^2/^1)063 = 0.45 ± 0.06 (solid line). 
We find that the radius of the I-band emission region 
in MG 04144-0534 is smaller than 0.70r?o with a statis- 



5 CONCLUSION 

Anomalous flux ratios between close pairs of images have 
been observed in a number of lensed quasars. SW02 showed 
that including a large smooth matter percentage in lens 
models can significantly increase the probability of these 
anomalous flux ratios occurring, if the source radius is small. 
We confirm this result, but also find that increasing the 
source radius destroys the broadening of the magnification 
distributions in lenses with a smooth matter comp onent. 
This result is also obtained bv lCongdon et ahl (|2007l ). for a 
different region of k — 7 parameter space. 

Using an observed I-band flux ratio of (^2/^1)065 ~ 
0.45 ± 0.06 for the le nsed quasar MG 0414-f-0534 
l|Schechter fc Moorel Il993l ). we place an upper limit 
(95%) on the size of the emission region of 2.62 x 
10^^ hjQ^^ {M /MqY^^ cm. The smooth matter percentage in 
the lens model was allowed to vary between 0% and 99%, 
and remains unconstrained. This limit assumes a constant 
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mass for the microlenses, and appears to be dominated by 
the source radius prior rather than the data. 

Our size hmit for MG 0414+0534 is consistent with 
limits on quasar continuum emission regions determined 
using other methods. For e xamp le, in the lensed quasar 
Q2237+0305, IWvithe et all l|2000l ) placed upper and lower 
limits (99%) on the R-b and emission region of 6x l O^^cm and 
2 X 10^'' cm respectively. I Wambsganss et al. ^ ('1990) placed an 
upper limit on the Q2237+0305 optical emission region of 
approximately 2 x lO^^cm. 

MG 0414+0534 is not the only lensed quasar displaying 
anomalous flux ratios in close images. PG 1115+080, SDSS 
J0924+0219, WFl J2026-45 36, and HS 0810+2 554 are four 
other well known examples (jPoolev et all 120061 ). An analy- 
sis similar to the one conducted here was undertaken for the 
broad emission line flux ratios in SDSS J0924+0219 using 
HST data (Kceton et al. 200ji). The authors concluded that 
acceptable models exist with the broad line region as large 
as ~ 2.3 X 10" (M/Mn)^/^ cm for a smooth matter percent- 
age of 80-85%. IPoolev et al.l (|2006l ) found X-ray flux ratio 
anomalies in PG 1115+080 to be a factor of 6 more extreme 
than their optical counterparts, and concluded that the op- 
tical continuum emission region is 10-100 times larger than 
expected from a thin accretion disk model. 

By extending observations of MG 0414+0534 across a 
number of filters, we would be able to conduct our anal- 
ysis for a series of different emission regions. Theories of 
quasar accretion suggest that different wavelengths are emit- 
ted from different regions in the source. Therefore with 
multi-band observations we could separate these regions in 
radius and directly probe the structure of the quasar central 
engine. 

In conclusion, we find that differential magnification in 
lensed quasar images shows that microlensing is important, 
but that it does not measure the smooth matter content 
of the lens. The size of the emission region and parity of 
the images are more significant factors than the smooth 
matter component. We place a 95% upper limit on the 
size of the I-band emission region in MG 0414+0534 of 
2.62 X 10^^h~g^'^{M/MQy^^cm. Unlike most previous anal- 
yses, this limit is independent of the unknown transverse 
velocity of the source. 
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